Detection of Tiny Mechanical Motion by Means of the Ratchet Effect 
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We propose a position detection scheme for a nanoelectromechanical resonator based on the 
ratchet effect. This scheme has an advantage of being a dc measurement. We consider a three- 
junction SQUID where a part of the superconducting loop can perform mechanical motion. The 
response of the ratchet to a dc current is sensitive to the position of the resonator and the effect 
can be further enhanced by biasing the SQUID with an ac current. We discuss the feasibility of the 
proposed scheme in existing experimental setups. 



I. INTRODUCTION 

In recent years the research on nanoelectromechani- 
cal systems (NEMS) developed as a broad and rapidly 
growing field 1,2 . The low mass and small size of NEMS 
make them excellent candidates for high-precision forced, 
mass^ and position detection 3 - » £ & 7 - & £, with a sensitiv- 
ity that is ultimately limited only by quantum me- 
chanics. Measurements on SiN— and GaAsii doubly- 
clamped beams coupled to single-electron transistors 
have achieved the sensitivity of four times above the stan- 
dard quantum limit. Moreover, NEMS have been re- 
cently envisaged as integrated elements into quantum in- 
formation circuits, such as on o 12 ' 13 i 14 or several^ qubits, 
opening up promising new prospectives also in the field 
of quantum information processing. 

A crucial issue for exploiting the remarkable features 
of NEMS is the detection of their position. To this pur- 
pose two different techniques are currently applied: Op- 
tical detection based on coupling of an oscillator to a 



microwave cavitj- 



16,17,18,19,20,21.22.2:! 



is very sensitive, al- 



though somewhat difficult to integrate in nanocircuits. 
On the other hand traditional electrical detection 2 ^ based 
on the oscillator back-action offers the advantage that 
excitation and detection occur in the same circuit. How- 
ever, this approach leads to lower sensitivity and does 
not allow an easy discrimination between electronic and 
mechanical excitations. 

Typically all existent detection schemes are based on 
high frequency measurements. In this article we suggest 
a new way to detect mechanical oscillations which avoids 
these difficulties and allows to detect the mechanical mo- 
tion with a dc measurement. Our idea is based on the so- 
called ratchet effect 2 ^. In very simple terms the ratchet 
effect can be explained by considering a particle moving 
in an external periodic potential and subject to a peri- 
odic (or random) force with zero average. If the potential 
breaks spatial inversion symmetry, a directed motion is 
possible even if the system is unbiased. By detecting the 
directed motion of a particle it is in general possible to 
gain information on the oscillating force. Here we discuss 
in details the case when the nanomechanical resonator is 
a part of a SQUID loop. It is well known that SQUIDs 
are currently used as ultra sensitive magnetic field detec- 
tors, and that they represent a solid-state implementa- 



tion of qubits based on Josephson effect 26 . A remarkable 
attention in the last few years has also been devoted to 
coupling SQUIDs to mechanical oscillators for detection 
or cooling ^ 13 ' 27 ! 28 ' 29 and, very recently, the first detec- 
tion of mechanical motion of a micromechanical resonator 
embedded in a dc SQUID was reported 3 ^. 

It has been shown 3 -^ that SQUIDs can behave as ratch- 
ets, where the superconducting phase difference at the 
junction 7 and the Josephson potential play the role of 
the coordinate and the ratchet potential, respectively. 
When such systems are biased by an ac current, a con- 
stant drift in the phase drop occurs under appropriate 
conditions, causing a finite dc voltage. 

In this work we propose to exploit the rectifying prop- 
erties of these SQUID ratchets in order to detect the 
motion of a nanomechanical resonator; indeed a mechan- 
ical oscillator embedded in a SQUID circuit threaded by 
a magnetic flux produces an ac signal across the SQUID 
that plays the role of a zero-average driving force, affect- 
ing the dc characteristic of this ratchet circuit. We find 
that the current-voltage curve is qualitatively changed 
by the mechanical oscillations with respect to the case 
without a moving part; in principle certain features can 
be extracted that can only be attributed to the presence 
of mechanical oscillations, such as Shapiro steps and a 
non-monotonic dip-and-peak structure in the symmetric 
part of the curve coming from them. This reduces the 
problem of the detection of motion of a suitable class of 
NEMS to an accurate measurement of the dc current- 
voltage characteristic of this device, avoiding the difficul- 
ties in high-frequency ac measurements arising in other 
detection schemes. 

The paper is organized as follows. In Sec. [TTJ we in- 
troduce the model and discuss the physically relevant 
scales. The results for the current-voltage characteristics 
are then presented; depending on physical parameters of 
the junctions employed, the adiabatic limit of a slowly os- 
cillating resonator may be of relevance and it is presented 
in Sec. IIIH However, higher coupling efficiency may re- 
quire slower time scales for the motion of the resonator, 
thus necessitating a more general discussion for arbitrary 
resonator frequency, which is done in Sec. lIVI We discuss 
the role of thermal fluctuations in Sec. [V] and we then 
proceed with a critical assessment of the experimental 
realization of this scheme in existing setups (Sec. IVT[) by 



2 



comparing the numerical values of the relevant physical 
parameters. Concluding remarks and a thorough descrip- 
tion of this detection technique are presented in Sec. lVIII 




5 10 15 20 25 

7 



Figure 1: (a)-(b) A model and the corresponding scheme of 
the ratchet circuit considered here. Three Josephson junc- 
tions are coupled to a mechanical resonator oscillating in the 
plane of the circuit, threaded by a magnetic flux, and biased 
by a current. The junctions have equal critical currents and 
capacitances but different shunting resistances, (c) The po- 
tential U (7) in the absence of a moving part (solid line) ; the 
average slope in the decreasing part is larger in modulus than 
the average slope of the increasing part, as can be seen by 
comparing with a cosine curve (dashed line). 



II. THE MODEL 

The circuit that we analyze is schematically depicted 
in Fig. [TJ It consists of a Josephson junction (JJ) and 
a mechanical resonator connected in parallel with two 
other JJs. The Josephson junctions are in the classical 
regime and they can be described with the RCSJ model, 
i.e. 



— = sxaji 

J-r 



2eR l I, 



■7* + 



fiC 
2e7 



■7» 



i = 1,2,3 



(1) 



In Eq.([IJ the current flowing through each junction is 
denoted by Ii, and 7^ is the related gauge-invariant 
phase difference. In Eq. (JTJ) C denotes the capaci- 
tance, and Ri the resistances, chosen as R\ = R and 
R2 = R3 = i?/2 (see Fig. [TJ. For simplicity we assumed 
that the three junctions are characterized by the same 
critical current I c . By considering overdamped junctions, 
u p iRC <C 1 {lo p i — ^f2eI c /TiC is their plasma frequency), 
the terms involving second-order derivatives of the phases 
can be neglected in Eq. ([JJ. One can then show that, if 



72 = 73 holds at a certain time, such condition is main- 
tained. We shall henceforth set 72 = 73 = 7/2, whereas 
the phase 71 is related to the other ones by the equation 



71 



71 = 0,1,. 



(2) 



ensuring a vanishing total phase difference along the loop. 
In the previous equation <E> denotes the total magnetic 
flux threading the circuit and $0 is the elementary quan- 
tum of flux. Neglecting the self-inductance of the loop, 
the flux is determined by the total area of the circuit, 
which in turn depends on the position of the resonator: 
<f> = B(A + lX(t)) (I is the effective length of the res- 
onator, X(t) is the position of its center of mass and A 
is the area of circuit when X(t) = 0). In the following 
we will assume that the motion of the resonator is fixed 
externally and equals X(t) = X cos(cut) (see Sec. IVI|) . 
Indicating by lb the current biasing the device, current 
conservation 2&(t) = h(t) + hit) leads to the following 
equation of motion for the phase 7, 



A dU 
^-97- 



I b (r) , irLBX Q 



$0 



fisin(Or) , (3) 



where r = uj*t is a dimensionless time variable, oj* — 
eRI c /h and fl = ui/cu*. We will consider external bias 
of the form ife(r) = Idc + lac cos(51 ac r), with I e i c denot- 
ing the dc component, and I ac and fl ac the amplitude 
and frequency of a monochromatic ac component. The 
periodic function 

U(7, f(r)) = -2 cos( 7 /2) - cos( 7 + + /(r)) (4) 

is the potential leading to the ratchet effect as dis- 
cussed by Zapata et aP 1 . In addition to the dimen- 
sionless external magnetic flux <f> = 2irAB /<&o, the po- 
tential U now includes also a time-dependent fluctuation 
part /(t) = (2irLBXo/$o) cos(17r) which arises from the 
oscillations of the mechanical resonator. The coupling 
between the SQUID and the mechanical resonator can 
be parametrized by the dimensionless parameter 



2ttLBXq 
$0 



(5) 



At finite temperatures the r.h.s. of Eq.Q will also include 
a noise term. 

Depending on the value of the external flux, the poten- 
tial term defined in Eq.(U) breaks 7 inversion symmetry 
and therefore leads to the ratchet effect, fn Ref. |3jJ this 
effect has been shown by biasing the SQUID with an os- 
cillating current leading to a finite dc voltage. Here we 
show that the ratchet effect allows one to detect the oscil- 
latory motion of the oscillator which enters the dynamics 
of the systems both in the potential through the fluctu- 
ating part of the flux and in the third term on the r.h.s 
of Eq.Q. We emphasize that, in order for the ratchet 
effect to appear, the system needs to be out of thermal 
equilibrium 2 ^. 
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Figure 2: (a) Dimensionless dc voltage drop vq — Vo/RI c as a 
function of the dimensionless bias current Idc/Ic, for different 
values of the flux. No oscillating mechanical part is present 
(a = 0). (b) The even part Vq of the same curves plotted 
for positive values of Idc only. Vq is zero for cj> = nir and 
v e Q {-<t>) = -v^). 



The dc voltage drop across the circuit can be expressed 
through the Josephson relation as 



h /&y\ ri c 

VQ = Te\-dt/ = — {dTl) 



(6) 



where the brackets (. . .) indicate time average (and possi- 
bly thermal average). A non- vanishing voltage drop can 
be obtained when the ratchet potential is asymmetric 
in 7, which occurs for values of the magnetic flux <fi ^= nir. 
In particular one finds that, for sufficiently large values 
of I ac , a finite dc voltage Vq ^ arises even in the pres- 
ence of a vanishing average drive (Idc = 0). Furthermore, 
for Idc t^ 1 the dc V — I curve is not and odd function 
of the current. It is therefore convenient to analyze its 
symmetric (even) part defined as 



Vo(Idc) + Vo(-Idc) 
2RL 



(7) 



Mechanical oscillations lead to two effects in the dy- 
namics of the SQUID. On one hand a term qualitatively 
similar to an ac bias arises, as shown by the last term in 
the r.h.s of Eq.Q. Secondly, a fluctuating term /(r) ^ 
appears in the potential, affecting the symmetry of the 
ratchet potential as a function of 7. In the presence of the 
resonator the analysis is thus more subtle. As we shall see 
in Sec. IIII1 the fluctuations of the potential tend to wash 
out the ratchet effect; although for some specific values 
of the parameters one can still obtain Vo(Idc = 0) 7^ 0. In 
most cases the ratchet effect survives only in some non- 
trivial features in the symmetric part of the V — I curve, 
which we will mainly focus on in the rest of the paper. 

Apart from some specific cases the equation of motion 
of the SQUID, ©, does not allow for an analytic solution. 
Most of our results are based on a numerical integration 
of the equation of motion. As we are interested in using 



the SQUID as a detector of the mechanical oscillations, 
it is appropriate to first briefly discuss the case where 
no oscillator is present, in order to highlight the differ- 
ence due to the mechanical motion. When the SQUID is 
biased by a dc current [h(T) = Idc] and there is no res- 
onator (a = 0) the voltage drop can straightforwardly be 
computed as Vb = 4w(RI c /2)/T, where An is the period 
of the potential (H|), and 



T = 



Idc/Ic - sin(7/2) - sin(7 + 0) 



(8) 



is the time required for the phase to span such an an- 
gle. The resulting V — I curve of the device is shown 
in Fig. [Ha). For an asymmetric potential ((f) 7^ nir) 
the curves are not odd in the current, and a non van- 
ishing even component arises [see Fig. C^b)]. This 
component will be shown to be affected by the pres- 
ence of the mechanical resonator. As the flux varies, the 
curve vo(Idc) shifts to the left and then to the right with a 
period of 2tt; this reflects in the amplitude and sign of Vq. 
Indeed, from Eq. © one can show that Vq(— (j>) — —Vq(4>) 
and v§0 = 0) = 0. 



III. SLOWLY OSCILLATING RESONATOR - 
THE ADIABATIC LIMIT 

We now discuss the role of the resonator. Before ad- 
dressing the full numerical solution of Eq. ([3]), we dis- 
cuss here the limiting situation in which O 1, i.e. the 
oscillator's frequency is much smaller than the plasma 
frequency of the SQUID. As we shall see, this regime is 
indeed suggested by the typical experimental conditions, 
where ft ~ 1CU 3 , indicating that the mechanical oscil- 
lations are much slower than the dynamics of the junc- 
tions. Moreover the same condition can be realized for 
the frequency of the ac component of the bias current. 
Under these conditions, the equation of motion can be 
solved within the adiabatic approximation. One can de- 
termine v (Idc) by evaluating, for a fixed position of the 
oscillator X(t) and value of the bias current Zb(r), the 
time T(t) required for the phase to span an angle 4-7r. 
Averaging over r one then obtains 



vo = 



\{^/T(r)) 



(9) 



(h. 



^ h(0)/lc ~ d 7 U( % f(6)) + f sin(0) 

Equation ^ will be analyzed in two different cases in 
which the bias current is constant or has also an addi- 
tional monochromatic ac component. 



A. dc bias current 

We start by analyzing the situation of a dc bias current, 
/b(r) = Idc- In Fig.[3Ja) and[3Jb) the even component (J7]) 
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Figure 3: Vo(Idc) curves for different values of the coupling 
parameter a. The two plots refer to different values of the flux: 
(a) <j> — 7r/2 and (b) (j> = 7r/20. While for small coupling one 
recovers the curves of Fig. ([2]), for high coupling the features 
are washed out. The effect is more pronounced for small flux 
bias. 



of the voltage is plotted as a function of Idc, for different 
values of the coupling strength a and flux <p. For small 
coupling a <C 1, we recover the exact results discussed in 
Sec|TT] [see Fig. |2Jb)] . On the other side, at high values of 
the coupling the main effect of the mechanical oscillations 
is to wash out the effects of the asymmetry of the ratchet 
potential. This is due to the fact that in this system the 
resonator oscillations also affect the ratchet potential via 
the term f(t). Notice that the disappearance of asym- 
metry occurs at lower values of a for a smaller flux </>, 
as one can see by comparing the dashed-dotted curves in 
Fig. [3£a) and^b). Thus the effects of the coupling to 
the mechanical resonator are more dramatic for <fi ~ nn. 
Note, however, that for the special values <f> = nn, where 
n is an integer, one still finds Vq — 0, as can be shown by a 
change of integration variables 7 — > — 7 and 9 — ► ir + 9 in 
Eq. ([9]) with Ifc(r) = Idc- We notice also that the experi- 
mentally accessible values of a are very small. Although 
in Sec lVII possible operative ways to enhance the coupling 
a will be proposed, these results lead to conclude that, in 
the presence of a purely dc bias Idc, it is difficult to gain 
information on a realistic mechanical resonator. A possi- 
ble way out is to inject current with a monochromatic ac 
component with frequency L) ac /2ir, small enough for the 
adiabatic approximation still to hold. This is discussed 
in the next subsection. 



B. ac monochromatic bias 

In this case some new effects that can be clearly at- 
tributed to the mechanical resonator indeed arise. One 
might naively expect that such effects simply originate 
from resonances at specific values of u> ac related to oj, 
yielding an enhancement of the effects due to the me- 
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Figure 4: Vo(Idc) curves for I ac — I c , Q ac — 2S1/3 and differ- 
ent values of the coupling parameter a (the same as in Fig. [3j . 
The two plots refer to different values of the flux bias: (a) 
<f> = tv/2 and (b) (f> = n/20. Again for small coupling the 
curves are indistinguishable from the uncoupled ones (a = 0) ; 
however for high couplings now a richer structure appear and 
its set-on occurs at a smaller coupling for <f> ~ mv. 



chanical oscillations. However, as we shall discuss below, 
the situation is more complex and the overall features of 
the current-voltage characteristics originate from differ- 
ent mechanisms. 

We assume I b (r) = I dc + he cos(f2 ac r) (Q ac = uj ac /uj*) 
and we restrict to the case VL ac = Og/p, with q and p in- 
tegers with no common divisor apart from 1. A straight- 
forward calculation allows to rewrite Eq. © as 



sin[7 + cj) + a cos(p9)] + sinfj^) 



(10) 



and to determine how the ac bias affects the V — I curves 
of Fig. 03 The result is shown in Fig. [4] for two different 
values of the flux. 

The even component Vq is suppressed for strong cou- 
pling, similarly to the dc bias case. However, a richer 
structure due to the interplay of the two components of 
the driving bias is observable. Another noteworthy differ- 
ence with respect to the dc case is that, at small values of 
the dimensionless flux, the suppression occurs at higher 
values of the parameter a, as one can observe by compar- 
ing the curves of Fig. [31(b) with Fig. EJb). This indicates 
that the presence of an ac current makes the even com- 
ponent of the characteristic curve much more robust to 
the coupling with the resonator. Such difference is par- 
ticularly striking at <f) = mr (n = 0,±1, . . .). For these 
particular values of flux Vq vanishes for a purely dc bias 
current, whereas a finite value of Vq is predicted for an ac 
current. Indeed the argument used to prove that Vq = 
if 4> — nix (a suitable change in the integration variables) 
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Figure 5: «o(^dc) for zero bias flux and different values of the 
ratio f2 ac /Sl. Here a = 1, while for a = the curves should 
be flat (v§ = 0). 



does not apply when I ac ^ 0, unless the frequency as- 
sumes specific values, namely those for which p and q are 
both odd integers. One can see the effects of the two 
combined oscillations in Fig. [5] Notice that these curves 
cannot be due to an ac bias alone, for if a = and I ac ^ 
the same argument of the previous section would apply, 
and the V — I curve at <j) = would be odd. 

These results suggest a possible operative detection 
protocol for the tiny motion of a mechanical resonator in- 
tegrated in a SQUID ratchet circuit: one may first mea- 
sure the V — I curve of the circuit biased by a purely 
dc current, for different values of magnetic flux <fi. This 
would allow to determine values of <j> such that the V — I 
curve is completely odd. Then, by adding an ac term to 
the bias current, without changing the magnetic field, a 
non-vanishing even component would arise in the V — I 
curve with varying the frequency of the ac bias. 



IV. ARBITRARY RESONATOR FREQUENCY 

In many situations the value of the coupling param- 
eter a is rather small (see Sec. I VI]) : for improving the 
coupling efficiency one then could in principle study de- 
vices with a higher value of the frequency f2, so that 
the coupling term ail cos(f2r)/2 can be increased. Un- 
der these circumstances the adiabatic approximation may 
not hold anymore and one should compute the current- 
voltage characteristic starting from the complete numer- 
ical solution of Eq. (|3|) . 

The results are shown in Fig. (Ha), where the V — I 
curves at strong coupling are plotted, in the absence of 
the ac bias current, at vanishing flux and for different 
values of the dimcnsionlcss frequency f2. One can clearly 
recognize the typical shape of the V — I curve for ac bi- 
ased Josephson junctions, characterized by Shapiro steps. 
However, while these features are usually attributed to 
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Figure 6: (a) Exact vo-Idc curves for different frequencies Q 
of the mechanical resonator. Here a — 1, I ac = and (j> — 
0. Even though no ac bias current is supplied, the system 
exhibits Shapiro steps, whose heights is related to Q. (b) The 
same but for a fixed frequency Q — 1 and different values of 
the coupling a. The width of Shapiro steps approaches zero 
for small values of a. 



an ac bias current, here they originate from the mechani- 
cal oscillations. Notice that the height of the steps varies 
with the proper frequency of the resonator, as it happens 
for Shapiro steps in isolated Josephson junctions. This 
suggests that in principle the frequency of the resonator 
can be directly read out by inspection of the height of 
the steps in the V — I curve of the system. The width of 
the steps, instead, is usually related to the amplitude a 
of the ac biasing signal. This is illustrated in Fig. HJb) , 
where one can see that, as a approaches zero, the steps 
become narrower and the V — I curve becomes smoother. 

Important differences emerge with respect to the adi- 
abatic limit. In particular the even component of the 
V — I curve does not vanish for <f> = nir and I ac = 0, 
and exhibits an interesting structure: ranges where the 
voltage is independent of Idc are separated by sequences 
of maxima and minima, as shown in Fig. [7] Interest- 
ingly, the amplitude characterizing these features in the 
even component of the V — I curve seems roughly in- 
dependent of the coupling parameter a. This represents 
a promising effect for the detection of realistic mechani- 
cal oscillators. However, as it is already the case for the 
whole curve vo(Idc) (Fig- IS]), the width of the different 
peaks is suppressed for small coupling. 

The effect of an ac bias is to enhance the amplitude 
of these features; indeed the ratio between the heights 
of the features in Fig. [5] and those in Fig. [7] is about 50. 
The numerical solution indicates that, in the absence of a 
mechanical oscillator (a = 0), the V—I curve is odd (vq = 
0), as in the adiabatic limit. We conclude that at <p = nil 
a detected signal has to be attributed to the mechanical 
resonator only. Even though the main difficulty at low 
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Figure 7: The even component vq of the voltage drop as a 
function of the dc bias current Idc at (f> — and £7 = 1 for 
different values of the coupling a (the curves are displaced 
vertically). The amplitude of the features is the same for 
all couplings, but their widths becomes smaller for smaller 
couplings. 

coupling still remains, these effects may be enhanced by 
applying an ac bias. 

V. FINITE TEMPERATURES 

In this section we discuss the effects of thermal fluc- 
tuations on the dc voltage drop of the system. To be 
consistent with the model adopted so far, we shall focus 
on the classical regime and account for finite tempera- 
ture by adding a stochastic term (white noise) in Eq. ([3]). 
These fluctuations are associated to two sources of dissi- 
pation: the finite resistance R/2 of the circuit, giving rise 
to electrical fluctuations in the bias current, and the finite 
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Idc/ h 

Figure 8: The even component vq of the voltage drop as a 
function of the dc bias current Idc at </> = and I ac ~ Ic, for 
different values of the coupling a. Here Q — 1 and Q. ac — 0.67. 



quality factor Q of the mechanical resonator, related to a 
random force producing mechanical fluctuations. Electri- 
cal fluctuations are taken into account by adding a term 

AI b (r)=I c] l^^ar) , (11) 

to the bias current /^(t), where £(t) is a white- 
noise process with zero average and correlation function 
(£( T )£( T ')) = <H T — T ')- Mechanical fluctuations may be 
introduced via a stochastic force 

F(t) = J q^-»?(*) ( 12 ) 

acting on the resonator. Here m is the mass of the me- 
chanical oscillator and rj(t) another white-noise process, 
statistically independent of £. In the previous sections 
we assumed the motion of the mechanical oscillator to 
be undamped (Q = oo), so that the mechanical noise 
should be neglected. However for realistic devices this 
treatment is still a good approximation if 

k B T <^Q^muj 2 xl (13) 

with xq denoting the amplitude of the oscillation. For 
high Q there is a regime in which temperature is low 
enough to fulfill condition (fl"3"]) . but high enough to ful- 
fill fceT S> tkJ, so that the behavior is neither quantum 
nor noisy. In the following we shall assume that the me- 
chanical noise can be neglected, and analyze Eq. ([3]) with 

y = an/2 sin(ftr) + ^ + ^ cos(^ qc t) + V&£(t) (14) 

where 9 = T/T and T = hI c /Aek B ^ 12K • I c /l/iA. As 
one may expect, for a certain value of the coupling pa- 
rameter a there is a value of above which the peculiar 
features of the characteristic curve are lost. For deter- 
mining the temperature scale at which this happens, one 
has to know To, which is set by the value of the critical 
current. Devices with large critical currents are in prin- 
ciple the best for this detection scheme, since they will 
be less affected by temperature. For a critical current as 
large as 1mA, the temperature scale is 10 4 K. As an exam- 
ple, consider Fig. [SI which is a detail of the u§(id c ) curve 
in Fig. [7] for a = 0.1, compared with its finite tempera- 
ture counterpart. Thermal fluctuations tend to wash out 
the peculiar features of the dc characteristic; the temper- 
ature scale at which this happens is O.lmK for a high-/ c 
device, corresponding to 9 w 10~ 8 . 

VI. EXPERIMENTAL PERSPECTIVES 

By now top-down fabrication techniques can produce 
single-cristal doubly-clamped beams whose typical di- 
mensions are 1/zm x 0.1 /jm x 0.1/xm and whose frequency 
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Figure 9: Some details of Fig. [8]for a — 0.1, I ac — 1, fi = 1 
and r2 ac = 0.67; the curves correspond to different tempera- 
tures: \/e = (solid line) and V& = 10" 4 (dashed line). 

ranges from tens of MHa^ to 1 GHz 33 ; using Silicon Ni- 
tride or Carbide, the density is about 3 x 10 3 kg/m 3 and 
hence typical masses are of the order of 10~ 17 kg. Lower 
values can be obtained using carbon nanotubes, which 
can also have higher oscillation frequencies due to their 
high bulk modulus^. The oscillation can be induced by 
driving the resonator either mechanically (with a piezo 
element applied on the substrate and driven by an ac sig- 
nal) or electrically (for instance by capacitively coupling 
the resonator with a gate electrode); a 1/^m-long beam 
can be easily excited up to oscillation amplitudes of lnm. 
To estimate the corresponding coupling parameter a one 
has to keep in mind that the magnetic field causing the 
coupling cannot be too large, otherwise the device loses 
superconductivity; a typical upper bound is B < 0.1T. 
With this values, one gets a w 0.1, the value used for the 
plots in Fig. [9] 

The intensity of the back-action on the beam can be 
estimated by the Lorentz force acting on it in virtue of 
the current flowing through it: Fi — LIB ~ 1/xm x 
1mA x 0.1T = 10~ 10 7V; this is to be compared with 
mcj 2 x Q ~ l(T 17 kg x 4tt 2 (1GHz) 2 x l(r 9 m = 4 x 1CT 7 A. 
(At such amplitudes the force may be no longer linear in 
the displacement and it would probably be larger than 
this estimate) . The back action can thus be neglected as 
a first approximation. 

As for the adiabatic approximation, the value of fi — 
huj/eRI c depends on the resistance of the junction. For 
tunnel junctions the low-temperature resistance is due 
only to quasi-particlc tunneling and increases exponen- 



tially with decreasing temperature; at temperatures close 
to the critical temperature the resistance can be esti- 
mated with the normal state resistance and I c R n can 
be as large as lmV, corresponding to lo* ~ 10 12 Hz. In 
such situations the adiabatic approximation could apply. 
However the resistance of a Josephson junction can vary 
greatly and thus the validity of the adiabatic approxima- 
tion depends on the device. 



VII. CONCLUSIONS 

A SQUID ratchet circuit can in principle be used to de- 
tect the motion of a nanomechanical resonator integrated 
in it by means of dc measurements only. The characteris- 
tic current-voltage curve of this circuit exhibits Shapiro 
steps whose heights can be used to determine the fre- 
quency of the mechanical oscillations. Once this fre- 
quency is known, either directly with this dc technique 
or with others, one can bias the device with an ac signal 
whose frequency should be chosen in order to improve 
the effects of the mechanical oscillations; we suggest to 
use an ac frequency of 2/3 the estimated frequency of 
mechanical oscillations. Now, by varying the magnetic 
field through the device, the symmetric part of the char- 
acteristic curve should evolve periodically (the period 
being set by the flux quantum); the amplitude of the 
curve will be maximum for frustrated values of the flux 
($ e = (n + l/2)$o) and minimum for unfrustrated values 
(<f> e = n$o)- If there were no mechanical oscillations in 
the circuit, then the symmetric part of the characteris- 
tic curve should be zero for unfrustrated flux bias, since 
the potential experienced by the superconducting phase 
would be symmetric; everything that can be seen in this 
case is due to the interplay between the ac driving field 
and the signal due to the mechanical oscillations of the in- 
tegrated nanoresonator. The shape of the characteristic 
curve can be compared with the finite temperature re- 
sults we illustrated in this work (see for instance Fig. [9]). 
We emphasize that no ac measurement is required for 
this detection scheme. 
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